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Abstract 

We present LVG and non-local radiative transfer calculations involving the rotational and hyperfine 
structure of the spectrum of N2H" 1 " with collisional rate coefficients recently derived by us. The goal 
of this study is to check the validity of the assumptions made to treat the hyperfine structure and to 
study the physical mechanisms leading to the observed hyperfine anomalies. 

We find that the usual hypothesis of identical excitation temperatures for all hyperfine components 
of the J=l-0 transition is not correct within the range of densities existing in cold dense cores, i.e., a 
few 10 4 < n(H 2 ) < a few 10 6 cm . This is due to different radiative trapping effects in the hyperfine 
components. Moreover, within this range of densities and considering the typical abundance of N2H" 1 ", 
the total opacity of rotational lines has to be derived taking into account the hyperfine structure. The 
error made when only considering the rotational energy structure can be as large as 100%. 

Using non-local models we find that, due to saturation, hyperfine anomalies appear as soon as the 
total opacity of the J=l-0 transition becomes larger than ~ 20. Radiative scattering in less dense 
regions enhance these anomalies, and particularly, induce a differential increase of the excitation 
temperatures of the hyperfine components. This process is more effective for the transitions with 
the highest opacities for which emerging intensities are also reduced by self-absorption effects. These 
effects are not as critical as in HCO + or HCN, but should be taken into account when interpreting 
the spatial extent of the N2H" 1 " emission in dark clouds. 

Subject headings: line: formation : profiles — molecular processes — radiative transfer — ISM: clouds 
: molecules : abundances 



1. INTRODUCTION 

N2H 4 " was one of the first molecular ions d etected in in- 
terstellar space llTnadHens fc TnrneHII 97fl). The J=l-0 
line of this species has been extensively observed toward 
cold dark clouds and protostellar cores to get some es- 
timates of th e phys ical conditions of the gas (see, e.g., 
iRergin et al] Iffl: iTafalla et al.| J200l: iHotzel et all 

u 



^Ofll^Belloche fc Andrei {2004\ : iCaselli et alJ \2£ 

These observations indicate that N2H + is a good tracer 
of the highest density regions of dark clouds. It seems 
that N2H + is less depleted onto dust grain surfaces than 
CO and other molecular species. This is probably re- 
lated to the fact that N 2 , the chemical mother species of 
N 2 H + , is more volatile and condensates at lower temper- 
atures than carbon monoxide. In addition, the complex 
hyperfine structure of N2H + increases the odds to have 
at least one optically thin hyperfine line component to 
probe the innermost regions of these clouds. Therefore, 
this species is in principle an interesting tool to study 
cold dark clouds. However, a drawback has been the lack 
of collisional rate coefficients between N 2 H + and molec- 
ular hydrogen (or helium). The observational data in- 
dicate some hyperfine intensity anomalies that could be 
due to selecti ve collisional processes or to radiative tr ans- 
fer effects (seelGonzale z- Alfonso & Cernicharol lll99 3) for 
the case the analogous case of HCN hyperfine intensities 
in dark clouds). 
For current research, it is clear that astronomers need 
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to know the state to state collisional rates of N2H+ with 
H2 and He. This will be even more necessary for ALMA 
due to the much higher angular resolution and higher 
sensitivity observations that it could provide of proto- 
stellar cores in several rotational transitions of N 2 H + (up 
to J =9-8). 

N2H" 1 " has als o been detec t ed in warm molecular clouds 
Thaddeus (1977)) where the lines are 



(e.g. tu rner 

broader and very strong. In these objects only the hy- 
perfine structure due to the external N-atom could be 
noticed as the splitting produced by the internal N-atom 
is lower than the intrinsic line width. Nevertheless, in 
order to correctly model the N2H + intensities emerging 
from these clouds, astronomers need a complete set of 
state to state collisional rate coefficients for high tem- 
peratures. 

A detailed study o n mol ecular ions excitation was car- 
ried out by lOreenl (1 1 Q TfJl and a set of collisional rate 
coefficients for N 2 H + colliding with He was provided. 
In that work, the rate coefficients were given for tran- 
sitions among N 2 H + rotational energy levels. Recently, 
Dani el et all l)2005l) extended the previous study by com- 
puting a new set of collisional rate coefficients for transi- 
tions among hyperfine energy levels and using a new po- 
tential energy surface. The range of kinetic temperatures 
is between 5-50 K and in a future paper, rate coefficients 
for temperatures up to 300 K will be provided. 

The paper is organized as follows: Section is devoted 
to the spectroscopy of N2H" 1 ". In sectionElwe present the 
mBeaw».MB«aJl)«QoW#m.a- Large Velocity Gradient (LVG) 
model to discuss the treatment of molecular hyperfine 
transitions. Comparisons with the cases of HCN and HC1 
are made. In section f|| we present results from non-local 
radiative transfer models applied to N2H 4 " for different 
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cloud structures. 

2. SPECTROSCOPY OF N 2 H+ 

The energy levels of N2H + are characterized by the 
quantum numbers J (rotational quantum number), F\ 
that results from the coupling of J with I\ {F\ = J + 1\, 
where I\=l corresponds to the nuclear spin of the outer 
nitrogen), and F (F = F\ + I2, where 12 = 1 for the in- 
ner nitrogen) We have m aintained the symbols used in 
IDaniel et all lj2004 I2005T) . except for J. The external 
nitrogen nucleus induces the largest splitting since its 
coupling constants are larger than those of the internal 
nucleus. Following iGordv fc Cookl (Il984fl . the hyperfine 
energy levels can be found by diagonalizing the molecu- 
lar Hamiltonian H mo i = BJ 2 — DJ + H coup u ng , where 
B and D are respectively the rotational and centrifugal 
distortion constants of the molecule and H coup ii ng the 
effective nuclear coupling Hamiltonian. The Einstein co- 
efficients (A u ^i) are given by the equation: 



Aj Fi F-*J'F;F' - V JF 1 F^J'F[F' X Tj^Sj^F^j/F'fa) 

with the line strengths given by: 

SjfiF-^J'FIF' = [FiF[FF'\ [J,yA I ^ ^ ^ 2 I (2) 

where [x] stands for (2x+l) and {.} is the Wigner-6j co- 
efficient. The values adopted for rotational and coupling: 
consta nts were provided by L. Dore (see IDaniel et alJ 
(2004)) an d where determined following the method de- 
scribed in iCa.selli et alJ ill99 -">l. We assume a dipole 
moment for N?H + of 3.4±0.2 D, as derived experimen- 
tally bv lHavenith et all l|1990|h It is in excellent agree- 
ment with the val ue of 3.37 D derived theoretically by 
iBotschwinal JTT)S 1 : . The resulting frequencies and line 
strengths for the J =1-0 transition are given in Table 
Figure |2] shows energy diagrams of the J=l and J=0 
hyperfine levels and indicates the line strengths of the 
hyperfine transitions. Figure shows the J=l-0 hyper- 
fine components together with their labeling. 

Due to hyperfine interactions, there are 9 distinct en- 
ergy levels for J >1, 7 for J=l and 3 for J=0. It is worth 
noting that the three energy levels in J=0 are indistin- 
guishable from a spectroscopic point of view, as their 
energy splitting is less than 10 -6 Hz. Thus, although 
there are actually 15 allowed hyperfine transitions con- 
necting J=l to J=0, there are only 7 resolved features, 
usually labeled as: 110-011, 112-012 (112-011, and 012), 
111-010 (111-010, 011, and 012), 122-011 (122-011 and 
012), 123-012, 121-011 (121-010, 011, and 012), and 101- 
012 (101-010, 011 and 012). Hereafter, we call hyperfine 
component each one of 15 transitions, and set of tran- 
sitions is used to refer to each group among the seven 
groups of blended components, using the labeling indi- 
cated above. 

3. LVG MODELS 

3.1. Effect of collisional rate coefficients 

To date, two sets of collisional rat es are a vailab le for 
the rotational structure of N2H" 1 ": IGreenI l)1975l) and 



IDaniel et alJ ((2005) 3 . The differences between state to 
state rate coefficients increase with both A J and J as dis- 
cussed in Dan iel et al . ( 2005). In Figure are compared 
opacities and excitation temperatures obtained using the 
two different sets in the same LVG code. There is agree- 
ment within a few percent, and the main differences are 
mainly due to changes in critical densities. At 10 K, the 
rate coefficients for AJ = 1 agree within 20% and the 
critical densities for the radiative lines differ by a sim- 
ilar factor in the opposite way. In particular, the rate 
coefficient for the transitio n J=l -0 is 22% higher using 
the results of IDaniel et all ll2005j) and the derived criti- 
cal density is 22% smaller. In other words, the recently 
derived rate coefficients make this transition thermalized 
at slightly lower densities. 

For densities below 10 4 cm -3 the excitation tempera- 
ture of the J=l-0 line would be very close to the cosmic 
background temperature due to the high dipole moment 
of N2H+ and the high critical densities (see Figure EJ- 
However, for column densities as low as 10 12 cm -2 (left 
panels of Figure^, the opacity of the J=l-0 and J=2-l 
lines would be large enough to produce significant ab- 
sorption effects on the radiation emerging from the in- 
ner and denser regions of the clouds. It is difficult to 
estimate from emission measurements the abundance of 
N2H" 1 " in these low density regions. However, by observ- 
ing the J=2-l and J=3-2 lines it could be possible to 
assess the effect of the external layers of the cloud on the 
emerging profiles of the low- J lines of N 2 H + , and thus 
to indirectly determine the spatial repartition in N 2 H + 
abundance. This effect will be discussed in detail below. 

3.2. LVG treatment of hyperfine transitions 

The LVG model used for the simulations presented 
in this work takes into account all possible collisional 
and radiative transitions among hyperfine levels. Nev- 
ertheless, it does not account for local overlap, and par- 
ticularly, the 15 transitions associated with J=l-0 line 
are treated independently. In order to check the ef- 
fect of the hyperfine structure, we performed calculations 
with and without it for three species: HC1, HCN and 
N 2 H + . Figure shows, for the three species, the ratio 
of the summed opacities of the hyperfine components, 
t(1 — 0)=T,Ti, to the opacity of the unsplitted rotational 
transition, tr(1 — 0). This ratio will be referred to as 
R t =t(1 - 0)/tb(1 - 0). The range of r(l - 0) explored 
is 0.1 to 100. This figure reveals differences between the 
two approaches. Under LTE conditions (large volume 
density), or in radiative equilibrium with the cosmic ra- 
diation background (low volume density), the two ap- 
proaches lead to the same estimate of the total opacity 
because under these conditions the hyperfine levels are 
populated according to the statistical weights : 

[F] 

UjFlF = Vhh\ nj (3) 

3 It was incorrectly stated, in lDaniel et~aTI i20CTl2005h that the 
potential energy surface (PES1 used bv IGreenI 119751 was calcu- 
la ted under the electron-gas approximation. Indeed, as explained 
in lGreenl 11975ft . the PES was determined using the self-consistent 
field approach which gives accurate results for interacting systems 
where the charge-induced dipole term is prominent. This is actu- 
ally the case for the N2H + - He system and explain the agreement 
between the two sets of collisional rate coefficients for the transi- 
tions among the first rotational levels. 
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where [J/1/2] is the total number of hyperfine sub- 
magnetic levels for a given rotational quantum number 
J. Such a population scheme of hyperfme energy levels 
occurs when lines are optically thin or in the domains of 
low and high volume densities, as Einstein (A u ^i) and 
collisional (C„_,i) rate coefficients reduce similarly after 
summation over initial and final hyperfme levels: 

[F}C JFlF ^j, FlF , = [Jhh]Cj^j, (4) 

FiF[FF' 

[F]A JFlF ^j, F , F , = [JhI 2 }Aj^j, (5) 

FiF[FF' 

If n(H2) ~ Aj^ji /Cj^ji , both collisional and radia- 
tive processes compete in populating the energy levels. 
Figure 01 shows that in this case the excitation processes 
are less effective to populate the first excited rotational 
levels if the hyperfme structure is considered. Thus, the 
total population of the fundamental level J=0 is underes- 
timated in a treatment that only includes the rotational 
molecular structure. It gives rise to an underestimation 
of the opacity of the J=l-0 and J=2-l lines, and to an 
overestimate of the opacity of the J=3-2 transition of 
N2H + and HCN. For HC1 there is always an overesti- 
mate of the opacity. 

When the hyperfme structure is taken into account, 
both Einstein and collisional rate coefficients have the 
same order of magnitude than the corresponding values 
for the rotational structure. Thus, the amount of radia- 
tive de-excitation, in the hyperfme description, is glob- 
ally larger for the rotational lines. Then, this makes R r 
increase with the total opacity of the lines. As opacity 
increases, the different Einstein coefficients of hyperfme 
lines lead to different excitation conditions for each hy- 
perfme component, i.e., different excitation temperatures 
(see Figured- In addition, in an inhomogeneous cloud, 
hyperfme lines may be excited in regions where other 
lines still have low excitation temperatures. This effect 
is analyzed in section Q] with the help of non-local radia- 
tive transfer models. 

HC1 behaves differently with respect to HCN and 
N 2 H + due to its large rotational constant (B ~ 313.0 
GHz). The first excited levels are hardly populated at 
low temperatures. We can roughly estimate the number 
of rotational levels reached by collisions to be ~ksT/hB. 
Thus, for HC1, most molecules are in the fundamental 
J=0 level for T K =10 K and R T ~ 1. Although some 
discrepancies do appear for R T in high- J transitions, the 
opacities of these lines are so low at Tk =10 K that the 
emerging intensity is negligible. At low temperatures, a 
hyperfme treatment of HC1 is thus not necessary. 

HCN and N 2 H + have similar rotational constants (44.3 
GHz and 46.6 GHz respectively). For these two species, 
the excited rotational levels are efficiently pumped by 
both collisions and radiation. Hence, the hyperfme struc- 
ture has to be considered to correctly derive column den- 
sities and abundances. For the J=l-0 transition of N 2 H + 
the error induced in the opacity estimate by neglecting 
the hyperfme structure varies from 20% (r ~ 1) to 100% 
(r - 100). 

3.3. J =1-0 Hyperfine Brightness Temperature Ratios 

Hereafter, we call T l ex the excitation temperature of 
the hyperfine component i. LVG calculations show that 



there are two domains in the estimate of T l ex for J=l-0. 
When opacities of individual components are low (n <C 
1), the excitation temperature is the same for all hyper- 
fine components. In this case, the opacities are propor- 
tional to line strengths s< : n = s«t(1 — 0). This be- 
havior is similar to what is expected in the LTE regime. 
For higher opacities, photon trapping induces different 
behaviors for the hyperfine components. 

In the LTE approximation, where T eK is assumed to be 
the same for all lines, the brightness temperature ratios 
vary monotonously with the rotational line opacity : 

T" 1 - e- T < 

— = 6 

Tl* 1 - e~ r i v ' 

Therefore, the brightness temperature ratios will 
change from the optically thin case, where they are given 
by the statistical weights 1 : 3 : 5 : 7, to unity in the 
optically thick case. Figure shows the hyperfine ra- 
tios, in LTE and non-LTE conditions, of the different 
hyperfine transitions with respect to the thickest one, 
i.e., JFiF=123-012. As expected for high densities, the 
LVG calculations are coincident with the result of the 
LTE approximation. The LVG calculations show that in 
a large domain of the (n(H2), N(N2H + )) plane, i.e. n(H2) 
< 10 6 cm" 3 and N(N 2 H+) > 10 12 cm" 2 /(km s" 1 pc" 1 ), 
the ratios are smaller than expected in the LTE approx- 
imation. Furthermore, ratios below the lower possible 
limit predicted under LTE conditions, i.e. 1/7, 3/7, 5/7, 
occur when hyperfine levels are not populated according 
to their statistical weights. Hence, fitting observational 
data using the hypothesis of identical excitation temper- 
atures for all J=l-0 hyperfine components will lead to an 
erroneous determination of the total opacity. From our 
calculations it seems that using this procedure underes- 
timates t(1-0). 

To discuss the behavior of the ratios, it is convenient 
to define an averaged excitation temperature, T ave , us- 
ing LVG results for individual rotational lines, and their 
hyperfine components, as: 

t "" = w, "( 1 + eSe)) (7) 

with T = hv/k B and r(l - 0) = £ n. 

Figure shows the LVG excitation temperature for all 
hyperfine components of the J=l-0 transition divided by 
T ave . This figure also shows the ratio of individual hyper- 
fine line opacities to the expected LTE opacities. Figure 
shows individual excitation temperatures and opacities 
for all hyperfine transitions. It is clear from these fig- 
ures that the effects are only important for large column 
densities, i.e., for large total opacities as discussed above. 
Non-LTE effects tend to reduce the spread of the opac- 
ities of the different lines transferring opacity from the 
thickest to the thinnest ones. For the hyperfine transi- 
tions associated to J=l-0, and from Figures El and we 
find the following trends : 

• for JFiF=123-012 : Tl x > T ave and n < s,-r(l-0) 

• for lines with initial quantum number F=2 : T l ex 
~ Tave and n ~ Sjr(l - 0) 

• for lines with initial quantum number F=l or F=0 

: T l ex < T ave and n > Si t{1 - 0) 
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We note that the variations of excitation temperatures 
and opacities for the different hyperflne components are 
anti-correlated, i.e., an increase of T l ex is accompanied by 
a decrease of Tj , and inversely (see Figure . This fact 
suggests that non-LTE effects mainly induce variations in 
the population of the JF\F=QYl level (the lower energy 
level in the reference hyperfine transition) compared to 
the LTE values. As for moderate and large opacities the 
brightness temperatures T B * , are less sensitive than T\ x 
to variations in Tj, the non-LTE effects tend to increase 
T B for the JFiF=123-012 line and induce a decrease of 
T B for lines with initial quantum number F=0 or F=l. 
Thus, for a given H2 density and N2H" 1 " abundance, ratios 
derived in the LVG approximation are smaller than those 
obtained under the LTE approximation (see Figure HJ. 

4. NON-LOCAL RADIATIVE TRANSFER FOR N 2 H+ 

The LVG calculations discussed in the previous sec- 
tion could be a reasonable approximation to the emerg- 
ing N2H" 1 " intensities from dark clouds when line opac- 
ities are low. However, as the line strengths of the hy- 
perfine components are different, the density structure 
of the cloud, radiation scattering and/or radiative cou- 
pling between different cloud regions could affect the 
population of N 2 H + energy levels. In order to check 
the validity of the LVG approximation and the differ- 
ent assumptions made earlier in this paper to interpret 
the N2H" 1 " J =1-0 transition, we carried out non-local 
calculations using the code developed and described by 
IGonzalez- Alfonso fc Cernicharo! i 1 9931) . For the opacity 
range considered in this work the code provides reliable 
and fast convergence. 

We assume a cloud at 160 pc with an angular diameter 
of 30" which corresponds to a radius of 3.6 10 16 cm (i.e. 
0.023 pc), and we consider three different sets of models. 
In the first one we consider a core of uniform density. 
In the second set the central core is surrounded by an 
envelope with size 3 and 6 times the core size. Finally, 
the last set corresponds to a collapsing cloud. In all 
models the turbulence velocity was varied from to 0.4 
km s" 1 by step of 0.1 km s _1 . Finally, we adopt a kinetic 
temperature of 10 K in all models. 

In all figures related to this section the ordinate scale 
is antenna temperature obtained by convolving the cloud 
brightness temperature with the beam of the 30-m IRAM 
radio telescope as follows: half power beam widths of 27", 
13.5" and 9", beam efficiencies of 0.76, 0.59 and 0.42, and 
error beams of 350", 220", and 160" for the J =1-0, 2- 
1 and 3-2 lines respectively. The error beam is always 
larger than the size of the dense regions and thus the en- 
ergy entering the telescope radiation pattern through the 
error beam just accounts for a few percent of the total 
intensity at most. Moreover, the contribution to J=l- 
from the extended envelopes surrounding the cores is 
fully taken into account by considering the convolution 
with the beam of the telescope (main beam and error 
lobe). The only plots for which the emerging intensity 
has not been convolved with the beam correspond to the 
two lowest right panels of Figure In these cases, we 
aim to show the influence of the high density region on 
the excitation of the molecules in the low density enve- 
lope and this would have been hidden by performing the 
convolution with the telescope beam pattern. 



4.1. Uniform density cores 

For this set of models the N2H 4 " abundance has been 
varied from 4 10 -10 to 6.4 10~ 9 , and the volume den- 
sity from 2.5 10 4 cm -3 to 4 10 5 cm -3 , by multiplicative 
steps of 4. Figure |H1 shows the emerging line profiles for 
the J =1-0, 2-1 and 3-2 transitions for all computed den- 
sities n(H 2 ) and abundances X(N 2 H + ). As expected, the 
brightness temperature increases with n(H 2 ) due to colli- 
sional excitation, leading to the thermalization of the the 
J=l-0 hyperfine components. We see that the brightness 
temperatures of the J=l-0, 2-1, and 3-2 transitions are 
more sensitive to a variation of n(H2) in the range 2.5 
10 4 cm -3 < n(H2) <4 10 5 cm~ 3 than to an increase of 
X(N 2 H + ). This can be understood by looking at figure 
in this range of n(H2) the excitation temperatures in- 
crease rapidly. Moreover, for J=l-0 total opacities below 
20, the excitation temperatures depend essentially on the 
volume density as radiative trapping effects are negligible 
and do not induce significant departure for the different 
T\ x from a single value. Consequently, in the optically 
thin case and for a given volume density, the line inten- 
sities are proportional to X(N2H + ) (raw corresponding 
to n(H2) = 2.5 10 4 cm -3 on FigurelHJl. Radiative trap- 
ping and line saturation only become important for large 
N 2 H + column densities. On Figure|Hl for the J=l-0 line, 
and on the panel corresponding to n(H2)=10 5 cm -3 , we 
see that increasing the abundance from 1.6 10~ 9 to 6.4 
10 -9 modify the ratio between the the JF Y F = 101-012 
and 121-011 sets of transitions. This change is due to ra- 
diative trapping and entails variations in the T\ x for the 
involved hyperfine components. Thus for large opacities, 
i.e. r(l-0) ~ 20, the excitation temperatures start to be 
sensitive to X(N 2 H+). 

One of the main results that can be derived from 
Figure El is that hyperfine intensity anomalies appear 
when increasing the N 2 H + abundance. This is sim- 
ilar to the effect found in HCN a nd discussed by 
IGonzalez- Alfonso fc Cernicharo! l) 1993th For total opaci- 
ties in the J =1-0 transition above 10, differential radia- 
tive trapping across the cloud starts to affect the rel- 
ative intensities of the hyperfine components. More- 
over, in an inhomogeneous cloud we could expect to 
be sensitive to different spatial extents for e a ch hy - 
perfine component (see [C ernichar o^ fc Guelinl (.1987); 
IGonzalez- Alfonso fc Cernicha.ro! l)1993D L This is even 
more striking when considering transitions with the same 
initial quantum number F, expected to have the same 
brightness temperature in the LTE limit. The most af- 
fected hyperfine set of components is the 121-011 set 
which indeed appears weaker than the 101-012 and 111- 
010 sets (see section 2). A low brightness te mperature 
for thi s set has already been observed (see e.g. iLee et all 
(1999): lCaselli et all IT9!)5 : :. We find the general trend 
Tb(IH-OIO) > Tb(101-012) > T B (121-011) and T B (U2- 
012) > Tb(122-011). This behavior originates from dif- 
ferential non-LTE effects in the hyperfine components of 
the J=0 F\=l levels (F=0,l,2). This is not surprising 
as the three transitions from J=l F=l (Fi=0,l,2) have 
different line strengths (see TableOJ. The 121-011 tran- 
sition mostly depends on the J F\F =010 level and is in- 
dependent of J FxF =012 while the 111-010 and 101-012 
sets mainly depend on JFxF=012. Thus, the 121-011 
set weakness compared to the other two must be related 
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to an underpopulation of the JFiF=012 level that pro- 
duces higher excitation temperatures for the hyperflne 
lines involving it. As shown in Figure |H1 the effect ap- 
pears for high abundances, i.e., large opacities, while for 
low abundances the three sets have the same intensities. 
Hence, it seems that it is not related to a collisional effect 
but to a radiative one. The sum of the line strengths for 
lines connecting the J=l to the J=0 level is larger for 
J-Fi-F=012 than for the others. Thus, for large opacities 
this level will be more affected by radiative trapping than 
the others. 

For the J=2-l and J=3-2 lines the effect of the abun- 
dance on the hyperflne intensities is similar. Large 
abundances lead to larger opacities, differential radia- 
tive trapping between the hyperflne components, and 
to anomalies in the observed line intensities. For these 
levels, line overlap is really important as the hyper- 
flne components are closer in frequency and should be 
taken into account in a detailed modeling. In collaps- 
ing clouds the velocity field could allow to have radiative 
coupling across the cloud between hyperflne components 
with different frequencies. This effect has been found 
to be extremely important in molecu les such as HCN 
ijGonzalez- Alfonso fc Cernichard 11993ft . In the models 
shown in Figure |H1 we only consider the effects of opac- 
ity for each line without taking into account population 
transfer between levels due to radiative connection be- 
tween them. Hence, the effects shown in this figure for 
large abundances are purely due to the line opacities. For 
the J=2-l and 3-2 lines, the self absorption induced by 
the external layers of the cloud is extremely important. 
When mapping clouds, the J=2-l and 3-2 transitions 
have little intensity outside the central position. How- 
ever, although the excitation temperature of their hyper- 
fine components could be low for offset positions, their 
line opacities are still large enough to produce significant 
effects on the emerging profiles toward the densest re- 
gions. For example, the model in Figure |H| with n(H2) 
= 2.5 10 4 cm -3 has low emission for all lines. However, 
the opacities are large enough for the strongest hyperflne 
components of the J=l-0 and J=2-l lines to produce the 
above mentioned effects. Figure El shows the opacities of 
the lines for n(H2) = 10 5 cm -3 which roughly agree with 
those found in the case n(H 2 )=2.5 10 4 cm -3 multiplied 
by a factor 4. We see that the opacity of J=2-l will 
be even larger than that of J =1-0, and of the order of 
0.75/7.5 for the strongest component and for the low- 
est/largest abundance of N 2 H + . As shown in Figured 
the excitation temperatures of the J=l-0 hyperflne com- 
ponents vary strongly in the range 10 4 < n(H 2 ) < 10 6 
cm -3 . The same will occur for the J=2-l and 3-2 hyper- 
flne components. Consequently, strong absorption will 
be produced in the low density layers of a cloud with a 
marked density structure. 

4.2. Core/Envelope Clouds 

We have finally explored more complex structures for 
the cloud to check the effects quoted above. Figure ITfil 
shows the results for two different cases. Left and cen- 
tral panels correspond to a simple core+envelope struc- 
ture. For the line profiles shown in the left panels 
(T K =W K, n(H 2 ) core = 4 10 5 cm" 3 , n(H 2 ) c „„=4 10 4 
cm" 3 , X(N 2 H+) core/em; =4 10" 10 , R core =3.6 10 16 cm, 
R erll , =2 R C ore), the opacity in the envelope for the J=l- 



line is very small. The emerging profile is practically 
coincident with the one arising from the core. However, 
the J=2-l line is optically thick in the envelope and im- 
portant absorption is caused to the emission from the 
core. As the opacity of the different hyperflne compo- 
nents depends on the line strength, weak lines are not 
affected (components at extreme negative and positive 
velocities), while the central components are much more 
affected. The J=3-2 lines are practically not affected by 
the envelope. 

A different situation occurs if the abundance of N 2 H + 
and/or the physical size of the envelope increase. The 
line profiles reported in the central panels of Figure Eil 
(T K =W K, n(H 2 ) core = 4 10 5 cm" 3 , n(H 2 ) em) =5 10 3 
cm" 3 , X(N 2 H+) core/em) =6.4 10" 9 , R core =3.6 10 16 cm, 
Re™ =6 Rcore) show that the strongest hyperflne compo- 
nents of the J=l-0 line are strongly affected by absorp- 
tion (the effect is similar to the one obtained in the mod- 
els shown in the previous section, see Figure|Hl where ab- 
sorption occurs in the most external core layers where the 
excitation temperature has decreased with respect to the 
center). For the J=2-l hyperflne lines the presence of the 
envelope affect dramatically the emerging intensities and 
the shape of the line. The J=3-2 line is less affected since 
the density in the envelope is not large enough to pump 
the J=2 and J=3 levels efficiently. The core also has an 
effect on the excitation temperatures of the N 2 H + lines 
in the envelope. The right panels of Figure llOl show the 
emerging profile from the envelope at a position located 
30" away from the center. The top right panel shows the 
comparison of the emission from the envelope alone and 
the emission obtained when the core is also present and 
excites radiatively the N 2 H + molecules of the envelope. 
The intensity of the J=l-0 line has increased by a factor 
~ 2. The other two right panels show the radial intensity 
distribution (not convolved with the telescope beam) for 
the J=l-0 and J=2-l lines (integrated intensity over all 
hyperflne components). The "heating effect" is clearly 
seen as a significant increase of the integrated line in- 
tensity even at distances of 3-4 times the radius of the 
core. These effects should be taken into account when 
interpreting the radial distribution of N 2 H + from obser- 
vations of the J=l-0 line (which is much easier to detect 
than the high- J lines in these low density regions). 

4.3. Collapsing Clouds 

Over the past decades, mm continuum observations 
and star count analysis ifWard-Thomnson et alJ 119941) 
have revealed that the density structure of cold dark 
clouds consists of an inner region of nearly uniform den- 
sity and of an outer region where the density decreases 
as r _p with p ~ 2.0-2.5. In this section, we use this type 
of density profile with a power law index p = 2 in or- 
der to test the influence of the velocity field on emerging 
spectra. Figure ITT1 shows the profiles of the J=l-0, 2-1 
and 3-2 lines for a 180" diameter cloud (at a distance D 
= 160 pc) with n = 4 10 5 cm~ 3 , r = 20" and a uni- 
form abundance X(N 2 H+) = 2 10~ 10 . The velocity fields 
compared are linear functions of r with slopes of ± 1.5 
km s _1 pc -1 (the velocity at the outer radius is ± 0.1 
km s~ 4 ) . 

In non-LTE conditions, the line excitation tempera- 
tures follow the density profile and, in the present case, 
decrease outward. Thus, in the case of a linearly de- 
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creasing velocity field, photons emitted at the center of 
the cloud appear blue-shifted for the molecules at greater 
radii which then absorb the red wing of the central emis- 
sion (the opposite occurs for a linearly increasing velocity 
field) . This effect is important if the lines are optically 
thick. Thus, in typical interstellar clouds conditions, self 
absorption affects the J =2-1 transition (see figure EH for 
the opacities) while the emerging profiles of the J =1-0 
and 3-2 lines are similar with and without velocity fields. 

On Figure 111! we see that the effect of the infalling 
material is to enhance the blue wing of the J=2-l central 
component which is self absorbed in the static model, 
the red wing of the central component being reduced by 
self-absorption. In the case of a velocity field linearly 
increasing with r, the opposite occurs and emission in 
the red wing is enhanced. Note that because the J =2- 
1 central component is composed of blended hyperfine 
lines with different line strengths (see figure ITT|l . changes 
in the line profile are not similar with respect to a change 
in the sign of the slope. 

Also, we stress that this effect is similar to what hap- 
pens to rotational lines with large opacity, except that, 
in the present case, the J=2-l central component corre- 
sponds to hyperfine blended lines of small intrinsic line 
width. Thus, characteristic features appear even if the 
velocity gradient is weak. 

4.4. Rate coefficients 

In star forming regions, the relatively large abundance 
of H2 makes this species the main colliding partner for 
other molecules. Estimated values of the collisional rate 
coefficients of molecular species in collision with para- 
H 2 (j = 0) can be obtained from the ones calculated with 
He as the collision partner. The underlying approxima- 
tion is to consider identical cross sections values for the 
two colliding systems, and then apply a scaling factor 
of 1.37 to the rate coefficients in order to correct for 
the associated different reduced masses. However, by 
analogy with the results obtained on the collision sys- 
tems involving HCO + , one could expect the cross sec- 
tions values with H 2 as a collision partner to be 2-3 times 
larger than those with He depending on t he se lected ro- 
tational transition. As discussed in section l3~Tl the main 
effect introduced by higher collisional rate coefficients is 
to lower the critical density of the rotational lines, which 
has direct consequences on the determination of both 
density and molecular abundance. The hyperfine rate 
coefficients are also changed t hrough variations o f the 
ratio of the different (see iDaniel et all l2004fi . In- 
deed, these variations modify the relative importance of 
hyperfine rate coefficients among a given rotational line. 
In order to assess the influence of such variations, we car- 
ried out calculations where two distinct approximations 
are employed to determinate the hyperfine rate coeffi- 
cients. The first consisted on rate coefficients propor- 
tional to the statistical weights of the final levels of the 
transitions a nd the second was b ased on the IOS approx- 
imation (see Dan iel et alJ I2005D . Figure El shows the 
emerging profiles resulting from the three sets of calcu- 
lations for a cloud at T K = 10 K, n(H 2 ) = 2 10 5 cm' 3 
and X(N 2 H+) = 5 10" 10 . It is readily seen that the rel- 
ative intensities of J =1-0 hyperfine lines are similar and 
that the largest difference (~ 15-20%) is encountered for 
the thickest JF 1 F = 123-012 line. The three sets give 



the same results for the thinnest JF%F = 110-011 line. 
Moreover, the intensity in the J=3-2 line associated with 
the IOS set of collisional rates is enhanced compared to 
the two other sets. This is related to the high propen- 
sity of the AF = AFi = A J transitions which lead to a 
more efficient pumping of the high- J levels. Note that 
it entails higher T ex for the J=2— 1 hyperfine lines which 
reduce the self-absorption feature in the main J =2-1 hy- 
perfine component. Finally, one would expect the colli- 
sional rate coefficients calculated with H 2 to affect simi- 
larly the intensities in the J =1-0 line, through variations 
of the ratio of the opacity tensors P& . Thus, as discussed 
above, the most important effect introduced by rate coef- 
ficients calculated with H 2 would be to scale all the rate 
coefficients and thus modify the density and abundance 
estimates. 

5. CONCLUSIONS 

We have used numerical methods in order to investi- 
gate the excitation properties of molecules with a hyper- 
fine energy structure, focusing specially on N 2 H + . Our 
conclusions are summarized as follows : 

1. The hyperfine structure must be taken into account 
in the radiative transfer calculations in order to 
derive the total amount of molecules present in a 
given rotational level. When the energy structure 
of the molecule is restricted to its rotational energy 
structure, high-J levels are more efficiently popu- 
lated and the opacities of the low- J transitions are 
underestimated. The error increases with the col- 
umn density of N 2 H + . Errors as large as a factor 2 
could be induced if the hyperfine structure of N 2 H + 
is neglected. 

2. For the typical temperature of dark clouds, i.e. 
T K =10 K, and n(H 2 ) < 10 6 cm -3 the ratio of hy- 
perfine brightness temperatures derived in the LVG 
approximation are always smaller than the ones 
predicted using the LTE approximation. Thus, this 
latter method is inadequate to assess densities or 
column densities from observed N 2 H + lines for the 
typical conditions prevailing in such clouds. 

3. The assumption of equal excitation temperatures 
for all hyperfine components belonging to the same 
rotational transition is valid for low opacities and 
fails down for high opacities due to radiative pro- 
cesses. For low opacities, we do not find any differ- 
ence in the excitation temperatures of the different 
hyperfine transitions induced by collisional rate co- 
efficients. Our calculations indicate that for tem- 
peratures in the range 5-50 K and for all volume 
densities, the excitation temperature of these lines 
will be identical if the total opacity of the J =1-0 
line is small. We stress that such collisional ex- 
citation effects, often found in the literature, may 
not be invoked to explain the intensity anomalies 
reported for N 2 H+. 

4. Non-local radiative transfer results show that the 
J =1-0 intensity anomalies arise while the opacity 
of this line increases. Moreover, photon scatter- 
ing in the low density envelopes affect the J =1-0 
hyperfine lines differently and tend to reduce the 
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emerging intensity of the thicker lines whereas the 
intensity of the thinnest lines is not modified by 
the envelope. Thus, in order to reproduce the in- 
tensities of the J =1-0 hyperfine lines, the only valid 
method is a non-local computation of the radiative 
transfer. 

5. For the typical conditions of cold dark clouds, the 
J=2-l line is a good probe of low velocity fields. 
Nevertheless, this transition is difficult to observe 
from ground based observatories due to the high 
opacity of the atmosphere at 186 GHz, although 
some work can be carried out in dry weather con- 
ditions. 
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Frequencies, Einstein coefficients (A„j) and line strengths (s u i) for the hyperfine components of the J=l-0 transition 
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Fig. 1. — Hyperfine transitions associated to J=l-0 in the optically thin case. Each resolved frequency is associated to multiple 
transitions among hyperfine levels (see text), and the standard labeling of the lines is indicated. 



10 





















































































































T 1 

T 1 

1 1 

o 
> 


0.556 


CO 
CO 
CO 

o 


CO 
CO 
CO 

o 


I s - 

1 — 1 

o 


o 

LO 


o 

LO 

o 


I s - 

T 1 

o 


CO 
CO 
CO 

o 


CD 
LO 
LO 

o 


CO 

cv 
o 

o 


o 

CO 
CO 

c\i 


o 


T — 1 

o 


1.250 




> 


V N 










V 














> 










< > 





101 

121 
123 
122 

111 

112 
110 



010 
011 
012 

Fig. 2. — Line strengths of the 15 hyperfine components in J=l-0 transition. The thickness of the lines indicate their relative weight 
compared to the others. Lines strengths are normalized in such a way that, summing over all initial J=l levels, gives the degeneracy of 
the final J—0 levels, i.e. 5 for JFiF=012, 3 for J F\F =011 and 1 for JFiF=010. Thus, the sum over all 15 transitions gives the total spin 
degeneracy. 
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Fig. 3. — Excitation temperatures (T ea! ) and opacities (r) obtained with two different sets of coll isional rate c oefficients and a LVG code, 
for a temperature of T = 10K. SET 1 (solid lines) refers to the latest computed rate coefficients IDaniel et alJI2005Tl and SET 2 (dashed 
lines) to the previously available one (Green 1975). Left and right panels correspond to N2H+ column densities of respectively 10 12 and 
10 13 cm -2 . A systematic velocity field of 1 km s _1 has been assumed for the cloud. 
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Fig. 4. — Ratio between the opacities of the J=l-0 (solid), J— 2-1 (dotted) and J=3-2 (dashed) rotational lines, R T , determined with 
and without hyperfine description for three molecules presenting hyperfine structure : N2H+ (left panels), HCN (central panels), and HC1 
(right panels). For each molecule the column densities are such that, at T=10K, opacities of the J=l-0 lines are r(l — 0) ~ 0.1, 1, 10 and 
100. The column densities are respectively : N(N 2 H+) ~ 2.2 10 11 , 2.3 10 12 , 2.5 10 13 and 2.0 10 14 cm~ 2 /( km s_1 PC -1 ); N(HCN) ~ 3.7 
10 11 , 3.0 10 12 , 2.5 10 13 , and 3.4 10 14 cm" 2 /(km s" 1 pc" 1 ); N(HC1) ~ 8.5 10 11 , 6.9 10 12 , 7.3 10 13 , and 7.7 10 14 cm" 2 /(km s" 1 pc" 1 ). 
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Fig. 5. — Ratio of the brightness temperatures obtained from LTE calculations (long-dashed lines) to those from LVG calculation 
(continuous lines), for a temperature of T=10K. The reference transition for both calculations is ( J-Fi_F)=123-012. The abscissa corresponds 
to the N2H+ column density (from 10 11 to 10 16 cm _2 /(km s _1 pc -1 )) and the ordinate is the volume density n(H2) (10 3 -10 8 cm -3 ). 
Under LTE different transitions with the same initial quantum number F have the same brightness temperature. Prom top to bottom 
the ratios of transitions with initial quantum number F = 0,1, 2 are shown. These are compared to LVG results for the ratios Ts(110- 



011)/Tb(123-012), T s (lll-010)/T s (123-012), and T S (112-012)/T B (123-012). 
thermalized conditions, dotted (LTE) and continuous (LVG) lines converge. 



Note that for high densities (n(H2) > 10 7 cm 6 ), i.e.. 
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Fig. 6. — Ratio of excitation temperatures and optical depths obtained in the LVG approximation to the corresponding quantities in 
the LTE limit. The temperature is T=10K and N 2 H+ column densities are 3.9 10 12 (left), 2.5 10 13 (center) and 2.0 10 14 cm _2 /(km s _1 
pc _1 ) (right column). For each set of transitions, the opacities Tj and the mean excitation temperature T% x are obtained by summing 
over lines with same frequencies and t(1-0) and T a „ e are obtained by summing over all J— 1-0 hyperfine lines (c.f. eq. |7J . Hyperfine 
set of transitions are referenced using the standard labeling (see Section |2j. We see that the variation of the excitation temperatures is 
anti-correlated with the variation of the opacities. 
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Fig. 7. — Excitation temperatures, opacities and brightness temperatures obtained in the LVG approximation, at T=10K, for the 15 
hyperfine components of the J =1-0 line. The columns, from left to right, respectively correspond to N(N2H+)=2.2 10 11 , 2.3 10 12 , 2.5 10 13 , 
and 2.0 10 14 cm _2 /(km s _1 pc -1 ). These column densities correspond to 7=1-0 total opacities of respectively r(l — 0) ~ 0.1, 1, 10, and 
25. 
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has an uniform volume density, a radius of 3.6 10 cm, a turbulent velocity of 0.1 km s , and a kinetic temperature of 10 K. In each box 
the line profiles are shown for N2H+ abundances of 4 10 -10 (dotted), 1.6 10 -9 (grey), and 6.4 10 -9 (solid). Prom top to bottom the three 
panels for each transition correspond to volume densities, n(H2), of 2.5 10 4 (bottom), 10 s , and 4 10 5 cm" 3 (top). 
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Fig. 9. — Opacities for the J=l-0 (top), J=2-l and J=3-2 (bottom) lines of N2H+. The cloud has a uniform volume density of 10 5 cm -3 , 
a radius of 3.6 10 16 cm, a turbulent velocity of 0.1 km s -1 , and a kinetic temperature of 10 K. N2H+ abundance is 4 10 -10 (dashed), 1.6 
10~ 9 (gray) and 6.4 10" 9 (black). 
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Fig. 10. — Emerging profiles for the J=l-0, J=2-l and J=3-2 lines of N2H+ for a core/envelope geometry (left and central panels). 
The core has a diameter of 30" and a density of 4 10 5 cm -3 . The envelope has a diameter of 90"/180" and a density of 5 10 4 /5 10 3 cm -3 
for the left and central panels respectively. The N2H + abundance is 4/64 10 -10 (left/central panels). The kinetic temperature is 10 K 
and the turbulent velocity is 0.1 km s _1 in the core and the envelope. Dotted lines correspond to transitions emerging from the core alone 
while solid lines correspond to those arising from the core+envelope system. The right panels show the emerging spectrum at an offset 
position of 30" (pure envelope emission corresponds to the gray curve). The effect of the radiative excitation in the envelope due to the 
photons arising from the core is clearly seen in the two lowest right panels. 
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Fig. 11. — Emerging profiles for a cloud with different velocity fields: static (black line), increasing linear function (gray) and decreasing 
linear function (dotted). The lower right panel shows the central lines of the J =2-1 line, with the positions of the individual hyperfine 
components indicated by lines proportional to their line strength. 
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Fig. 12.— ,7=1-0, 2-1 and 3-2 emerging profiles (left column) for a cloud at T K = 10K, n(H 2 ) = 2 10 5 cm" 3 and X(N 2 H+) = 5 
10~ 10 . Two approximations for the hyperfine collisional rate coefficients (see text) are compared to the set of reference (black solid line). 
The first consist of rates proportional to statistical weights (gray solid lines) and the second is based on the IOS approximation (dashed 
lines). The right column corresponds to a blow-up of the J=l-0 hyperfine lines. 



